Manifestations of electron interactions in the photogalvanic effect in chiral nanotubes 
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Carbon nanotubes provide one of the most accessible experimental realizations of one dimensional 
electron systems. In the experimentally relevant regime of low doping the Luttinger liquid formed by 
electrons may be approximated by a Wigner crystal. The crystal-like electronic order suggests that 
nanotubes exhibit effects similar to the Mossbauer effect where the momentum of an emitted photon 
is absorbed by the whole crystal. We show that the circular photovoltaic effect in chiral nanotubes 
is of the same nature. We obtain the frequency dependence of the photovoltage and characterize its 
singularities in a broad frequency range where the electron correlations are essential. Our predictions 
provide a basis for using the photogalvanic effect as a new experimental probe of electron correlations 
in nanotubes. 
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I. INTRODUCTION 

Carbon nanotubes (CN) can be viewed as long cylin- 
ders made from a graphene sheet [iHsl. Depending on 
the way this sheet is rolled up, the cylinder may have 
chiral (i.e., helical) structure. Chirality implies the exis- 
tence of the circular photogalvanic effect (CPGE), where 
a circularly polarized electromagnetic wave induces a dc 
current. The magnitude of the photo-induced voltage is 
finite even when the momentum of the photon is negligi- 
ble. This effect was predicted in Refs. and since then 
was investigated both experimentally and theoretically - 
see, for example, Refs. [6|-|9| and references therein. 

Usually the circular photogalvanic effect in non- 
centrosymmetric media arises due to a spin-orbit inter- 
action. In carbon nanotubes the spin-orbit interaction is 
weak and the circular photogalvanic effect arises because 
the electron motion along the tube axis and circumfer- 
ence is strongly mixed in chiral nanotubes. The single 
particle theory of chirality related optical effects in car- 
bon nanotubes has been developed in Ref. ■ This the- 
ory holds for non-degenerate electrons in semiconductor 
nanotubes or for metallic ones at high enough frequencies 
of the electromagnetic field. 

However, at low temperatures the single particle de- 
scription of one-dimensional interacting electron sys- 
tems, such as metallic CN, is invalid. Being a quasi- 
one-dimensional system, CN exhibits distinctive features 
which result from the strong correlations of the electrons 
in the system ll|, for instance, the Fermi edge singu- 
larity which is manifested by a characteristic power-law 
singularity of the tunneling density of states. The pho- 
tovoltaic effect in this regime was not investigated, and 
will be studied here. 

The paper is organized as follows. In Sec. II we re- 
view the non-interacting theory of the CPGE in CN. In 
Sec. Ill we study the problem with electron-electron in- 
teractions. Next, in Sec. IV, we discuss the limit where 






"2 















(b) 



FIG. 1: (a) The honeycomb lattice of graphene and its prim- 
itive lattice vectors a and b. (b) The Brillouin zone of 
graphene. K and K' denote the Dirac points of the energy 
spectrum. The straight lines cut the zone at the discrete set 
of allowed transverse momenta. This set is determined by 
the way a graphene sheet is rolled into a tube, such that the 
wave function is invariant with respect to the translation vec- 
tor L = na + mb. The example shown here corresponds to 
(n, m) = (5, 12). The lowest energy excitation near the K and 
K' points of the Brillouin zone are associated with a change 
of the angular momentum by ±h, thus a circularly polarized 
light directed along the tube excites electrons only near one 
of the Dirac points. 



interactions are strong enough so that the electronic cor- 
relations are governed by Wigner crystal order. We draw 
parallels between the CPGE and the Mossbauer effect, 
and show that the photovoltage, as a function of the fre- 
quency of the electromagnetic wave, exhibits a series of 
distinct singularities. Our results are summarized in Sec- 
tion V. 



II. SINGLE PARTICLE PICTURE 

Let us begin with a brief review of the noninteract- 
ing theory of the CPGE in CN's. Fig. 1(a) presents 
the honeycomb lattice structure of a graphene sheet, 
and its two primitive lattice vectors a = a(0, 1) and 
b ~ a(— 1, -\/3)/2, where a = 2.49A is the lattice con- 
stant. A single- walled CN rolled up from such a sheet is 
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defined by two integers n and m which specify the trans- 
lation vector L = na + mb that wraps around the cylin- 
der. The requirement that the wave function be invari- 
ant with respect to a translation by L, ip{x. + 'L) — V'(x), 
imposes a quantization condition on the component of 
the wave vector in the transverse direction k±\L\ = 27rj, 
where j is an integer. Fig. 1 (b) shows the Brillouin zone 
of graphene. The corners of this zone are Dirac points 
where the spectrum is degenerate and the two inequiva- 
lent points are denoted by K and K' . The straight lines 
cutting through the zone represent the aforementioned 
quantization condition. Each line, i.e, each value of j, 
defines a subband of the nanotube spectrum which cor- 
responds to a definite angular momentum of the electron 
around the cylinder (see, for example, Refs. lol. 12. 13|). 

A characteristic feature of chiral nanotubes is that the 
minima of subbands with different values of angular mo- 
mentum are shifted with respect to each other by a mo- 
mentum 5p, as illustrated in Fig. 2. To be concrete we 
consider the situation where the lowest conduction sub- 
band with angular momentum Tij is partially occupied 
(say due to doping by an an external gate), and the an- 
gular momentum of the first unoccupied band of the K 
valley is ?i(j + l). Now let us assume that a circularly po- 
larized electromagnetic wave is radiated in the direction 
of the tube. Absorption of a photon will be accompanied 
by a unit change of the angular momentum, which at low 
enough frequencies might take place only in one valley, 
say K, since the corresponding transition in the opposite 
(K') valley is forbidden by conservation of angular mo- 
mentum; see Fig. 1(b). A finite value of the photovoltage 
emerges from the asymmetry of velocities and light ab- 
sorption probabilities for right- and left- moving electrons 
in the K valley. The single particle absorption threshold 
ojI^ is determined by the electron band structure and is 
illustrated in Fig. 2. Because of the constant density 
of states at the Fermi level, the photo-voltage V{uj) as a 
function of the frequency uj experiences a jump from zero 
to a finite value at the threshold frequency uj = uj^p. 



III. INTERACTING PARTICLES PICTURE 

Generally, the frequency Wc of the many-body absorp- 
tion threshold is different from the single particle thresh- 
old uj^P. At sufficiently large values of 5p the absorption 
threshold can be significantly lower. It corresponds to an 
indirect transition in which one electron is transferred to 
a state near the bottom of the upper subband. The to- 
tal momentum of the system does not change as a result 
of photoabsorption. Therefore a momentum —Sp is im- 
parted to the electrons in the lower (partially occupied) 
subband. Since the velocity of the excited electron in 
the upper band is small, only the electrons in the lower 
band contribute to the photocurrent. The photovoltage 
V{uj) necessary to nullify the current is determined by 
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FIG. 2: Subbands of chiral nanotubes demonstrating the 
asymmetric band structure near each valley. The shaded re- 
gion in the lower box represents the occupied states within 
the single particle picture. 



the momentum flux required to compensate for the fiux 
of momentum to the lower band electrons due to pho- 
ton absorption. Denoting the photon absorption rate by 
W{ijj), we obtain for the photo-induced voltage 



V{u:) = ^VF(w), 
ne 



(1) 



where n is the electron density in the lower band, and 
e is the electron charge. The precise criteria of validity 
of Eq. (HI) depend on the value of and will be dis- 
cussed below. Here = U/Ekin is the ratio between 
the characteristic potential, U ~ e^n, and the kinetic, 
Ekin ^ (hn)'^ /m (m being the electron mass) energies of 
the electrons. 

Another manifestation of the many body effects is 
that near the threshold, |w — <C cjc, the photovolt- 
age exhibits a power law singularity V{uj) ~ (w — uJc)^° 
characteristic of one dimensional (ID) systems. We will 
show that rjo > 0, for Sp ^ hnrl^^, and 770 < when 
6p <^ hnrl^'^. Thus interactions suppress the photogal- 
vanic effect near the threshold in the first case and en- 
hance it in the second. We shall also discuss additional 
singularities in W{uj) which occur above the threshold at 
intervals of the Debye frequency of the Wigner crystal 
and are associated with the generation of hard plasmons. 

In principle, recombination processes create an addi- 
tional contribution to the electron distribution function 
in the lower band which is asymmetric in p and which 
affects the photovoltage V{u!). In this case however, the 
singularity of V(uj) at the absorption threshold remains 
intact. Recombination with a large energy transfer usu- 
ally takes place through a multiphonon emission via deep 
impurities which are associated with a large momentum 
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transfer to individual phonons. Therefore the momen- 
tum asymmetry of the distribution function in the lower 
band associated with such processes is small and can be 
neglected. 

The significance of many-body effects for the absorp- 
tion and the photovoltaic effect can be appreciated even 
in the framework of perturbation theory using simple 
kinematic considerations. Let us assume that Sp is larger 
than the characteristic momenta of electrons in the lower 
band. Consider an absorption process in which one elec- 
tron is transferred to a state near the bottom of the 
upper band and N electrons in the lower band receive 
a recoil momentum of order p* « —Sp/N and energy 
e* ~ {p*)^/2m. The energy of such a photoexcitation 
is TiutN ~ A -I- {Sp}'^ /{2mN), where A is the energy dif- 
ference between the bottoms of the bands (Fig. 1). It 
decreases with increasing N and may be significantly 
lower than the single particle threshold lu^p ~ i^N=i- If 
e* ^ ma.x[U, Ep], regardless of the details of the elec- 
tron wave function in the ground state, the excited elec- 
trons may be treated as free particles and the absorp- 
tion rate W{u!) can be calculated using the iV-th order 
of the perturbation theory in a way similar to that of 
Refs. [13, llBl- The increase of N with decreasing fre- 
quency can be considered as a precursor of the Lut- 
tinger liquid regime. Practically, however, the condition 
((5p)^/(2m) ^ ma.x[U, Ep] is difficult to realize in car- 
bon nanotubes, and therefore there is no parametrically 
big interval of frequencies where the perturbation theory 
works. For this reason, below, we focus our attention on 
the case where frequency is close to the threshold. 



IV. THE WIGNER CRYSTAL LIMIT 

Consider the experimentally relevant regime of low 
electronic densities n, where the ratio of the interac- 
tion energy to the kinetic energy of the electrons is large, 
rs ^ I. In this case exchange processes associated with 
the tunneling of electrons are exponentially suppressed 
in Ts, and will therefore be neglected in what follows. In 
this approximation electrons may be labeled by the site 
number i/ in the Wigner crystal lattice. 

The wave functions of electrons localized near a given 
site are superpositions of Bloch functions either in the 
valley K or K' . The two possibilities are realized with 
equal probability. In a broad frequency range electrons 
excited by light into the upper subband remain localized 
at their original sites of the Wigner crystal. For a cir- 
cularly polarized light only electrons from valley K can 
participate in the absorption. Thus the wave functions of 
excited electrons are superpositions of Bloch wave func- 
tions of the upper band in the K valley. Since they are 
localized near the Wigner crystal sites, their average mo- 
mentum must be 6p; see Fig. 1. Conservation of mo- 
mentum dictates that the opposite momentum —Sp must 



be transferred into the collective motion of the Wigner 
crystal, which leads to Eq. ([1]). A similar transfer of 
the recoil momentum from a given atom to the whole 
crystal occurs during the emission of 7 rays. One of the 
striking manifestations of the transfer of the recoil mo- 
mentum to the collective motion in three dimensional 
crystals is the Mossbauer effect [l^ in which no phonons 
are excited during a 7-ray emission. The probability 
for such a transition is given by the Debye- Waller fac- 
tor exp[— where Sp is the recoil momen- 



tum and 



is the variance of the displacement of the 



nucleus from its equilibrium position. In a one dimen- 
sional system the latter diverges logarithmically with the 
system size and the Mossbauer effect is impossible. How- 
ever light absorption with the emission of collective exci- 
tations (plasmons) remains possible. The probability of 
plasmon emission determines the absorption rate above 
the threshold, u > Uc- Though these plasmons carry the 
excess energy, their average momentum is zero. Thus 
each photon absorption provides a momentum Sp which 
is transferred to the Wigner crystal of the electrons in 
the occupied subband. Since the average velocity of the 
excited electron in the upper band is zero, we again arrive 
at the expression for V{u!) given by Eq. ([T]). 

Near the absorption threshold the probability of inter- 
band transitions exhibits a power law singularity, which 
can be described phenomenologically in terms of a motion 
of a mobile impurity in a Luttinger liquid 17H23|. The 
strength of electron-electron interactions in the Wigner 
crystal regime enables us to determine the light absorp- 
tion probability not only near the absorption threshold 
but in a much wider frequency range. 

Although many of our conclusions have general charac- 
ter, we will consider the case {Sp)'^/2m ^ ne^ where the 
Wigner crystal may be viewed as a weakly anharmonic 
chain 24]. We model our system by the Hamiltonian 
H = Hq + Hep, where the unperturbed Hamiltonian is 



A- 



2m 





2m 



(2) 



Here p^ and a;^ denote the momentum and position of the 
I'-th electron, A is the energy shift between the bands, 
and V{x) is the interaction potential. For simplicity the 
electron mass is assumed to be the same in both sub- 
bands. The matrix in Eq. ([2|) acts on the subband index. 
In the rotating wave approximation the interaction with 
the electromagnetic field is described by 



h.c). 



(3) 



Here a is the photon annihilation operator, cr+ is the 
raising operator of the v-th. electron from the lower to 
the upper subband, and D is the dipole matrix element 
of the transition. The electron positions and momenta in 
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a Wigner crystal regime can be expressed in terms of the 
plasmon annihilation and creation operators, bq and b^g, 



— v/n + , 



2mNuj, 



(6, + 6%)e''^ (4a) 



(4b) 



where n and N are the electron density and the number 
of sites, respectively, in the Wigner crystal, and LOq is 
the frequency of the a phonon with (dimensionless) wave 
number q which stratifies the dispersion relation: 



(PV{x) 



[1 — cos{qv)\ 



(5) 



x—vj n 



Denoting by c and the annihilation and creation of 
photoexcitations, respectively, on site v — 0, one can 
write the Hamiltonian in the form 



Ho = rl- 

q \ 
+ Ac'c, 



iSp c^c 




iSp c 
^2mNhijjq 



(6) 



where we have subtracted the zero point energy of the 
plasmons. The absorption of photons happens indepen- 
dently at different lattice sites; thus, within the dipole 
approximation, it is given by the Kubo formula. 



W{uj)-- 



Re I dte^"*([c(i),ct(0)]), (7) 



where (• • • ) denotes thermal averaging. To decouple the 
interaction between the plasmons and the excited particle 
we perform a unitary transformation of the operators 



U = exp j i6p c'c'y^^ 



\ + hlq 



(8) 



^ ^ yj2mNTiujq J 

In the transformed basis the Hamiltonian takes the form 
U^HqJJ = ^ hujab\bq + hojjc, (9) 



while the creation and annihilation operators are 



U'^S ^ ce^^P""'''^, and tj'^c'^U 



ce 



-iSpxQ /h 



(10) 



Thus in the limit of zero temperature Eq. ([T]) reduces to 

{6p?F{t) 



/(it exp I i{Lu — Is. / K)t 
-oo \ 



ti'^ J-oo ' V ' ' ' 2mhuJD 
where wd = is the Debye frequency, and 



ai) 



(12) 




FIG. 3: The absorption spectrum obtained by numerical in- 
tegration of (lllf) for the case of short range interactions and 
a momentum shift Sp, for which ryo = —0.4 



The singular behavior of the absorption rate is deter- 
mined by the long time asymptotic behavior of F{t). In 
this limit the integral over q may be evaluated by the 
stationary phase approximation. For a typical disper- 
sion, the stationary phase is at g = tt and the dispersion, 
near this point, may be approximated as 



(13) 



with Lo'^ = cP'u!q/dq'^\q^T^. But there are also contribu- 
tions from the end points of the integration interval, at 
q — Q and q — 2n. Near these points, the dispersion may 
be linearized. In particular, near q — 0, 



(14) 



where 



= duiq/dq\q^o, and similarly near q — 27r, 
ujq{27t — q). Evaluating the integral ([T2|) using these 



approximations we obtain 



Fit) 



2uj 



D 



In(iajQt) -I- 7 ■ 



exp(— ia;£)i) 



(15) 



where 7 is a constant of order unity which depends on 
the precise form of ujq. The first term of this asymp- 
totic formula accounts for the soft plasmons generated by 
the excitation, while the last term is associated with the 
excitation of hard plasmons with energy near the plas- 
mon Debye frequency uj-q. Substituting this expression 
in (jlip and expanding the hard plasmons contribution 
in a power series we find that the absorption exhibits 
singularities at frequencies ujj = A/?i H- jwc which are 
associated with the generation of j hard plasmons. Thus 
for — ^ ujd we get 



where 



Sp-" 



TTmhujQ 



2 



(16) 



(17) 
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and Sj{uj) = aj + Pj9{uj) is a step function between two 
values which depend on j. 

The existence of higher frequency singularities, j > 0, 
Eq. (|17l) . is a result of the harmonic approximation, 
which is valid as — > cxd. At finite values of rg these 
singularities are smeared by anharmonic interactions be- 
tween plasmons. The power of the higher frequency sin- 
gularities, r]j, increases with j. Therefore only the first 
few of them are significant. An example of this behavior 
in the case where ujq = wd| sin((3'/2)| and the value of 
5p is chosen such that 779 = —0.4 is depicted in Fig. 3. 
The power law singularity at the absorption threshold 
(j = 0) is not broadened by anharmonic interactions. Its 
exponent is 770 « - 1. 



Depending on the momentum shift between the sub- 
bands, the photovoltage diverges or vanishes near the 
threshold. Finally we remark that the mass difference be- 
tween the subbands, and the electron-phonon interaction 
do not change our qualitative picture. The mass differ- 
ence will slightly increase the singularity powers rjj due 
to the sudden change of the effective interaction strength 
between the excited electron and the rest of electrons, 
which leads to a contribution to r/o associated with the 
orthogonality catastrophe. 

Let us turn now to the case ^ 1 where interactions 
are weak and electrons can move freely along the sys- 
tem. In this limit the threshold frequency approaches 
the single particle value uj^^ and usually Sp <^ pp 
where pp is the Fermi momentum of the occupied band. 
The photoexcitation induces a direct transition which 
does not change the electron momentum. Therefore the 
Mossbauer mechanism described above does not apply. 
However, if the energy relaxation time of an excited elec- 
tron in the upper subband is shorter than the momen- 
tum relaxation time of the electrons in the bottom band, 
then each photoexcitation results in an effective momen- 
tum transfer of —pp to the electrons in the lower band. 
Thus the photovoltage is given by Eq. ([T]) in which Sp 
should be replaced hy pp. The absorption rate has a 
power law singularity exponent at the absorption thresh- 
old, W{u!) ^ {uj — uj^P)^" with the exponent 770 ~ r^. The 
above threshold singularities shown in Fig. 3 are expected 
to smear out. 



V. SUMMARY 

In this paper we have shown that the inclusion of in- 
teractions into the CPGE problem in CN changes the 
frequency dependence of the photovoltage V{uj) dramat- 
ically. Instead of a jump at the threshold frequency, the 
interaction produce a singularity with a power exponent 
which depends on the relation between the momentum 
shift between subbands and the strength of the interac- 



tion between neighboring electrons, similar to the classi- 
cal Mossbauer effect. We also show that V{uj) exhibits 
additional singularities at higher frequencies. These re- 
sults provide a new experimental probe by which the na- 
ture of electronic correlations in carbon nanotubes can 
be examined. 
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